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Abstract 

B 

We theoretically investigate surface modifications on alkaline earth halides 

d 

due to highly charged ion impact, focusing on recent experimental evidence 

for both etch pit and nano-hillock formation on CaF2 pp. We discuss media- 
CD 

i Q | nisms for converting the projectile potential and kinetic energies into thermal 

energy capable of changing the surface structure. A proof-of-principle clas- 

> 

V.Q sical molecular dynamics simulation suggests the existence of two thresholds 

which we associate with etch pit and nano-hillock formation in qualitative 

in 

agreement with experiment. 
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•i-h 1. Introduction 
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Investigations of the interactions of slow [v < fBohr), highly charged ions 
(HCI) with surfaces were originally aimed at gaining an improved under- 
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standing of the dynamics governing the neutralization and eventual dissipa- 
tion of the large projectile potential energy in a nano scale region of the target 
on an ultrashort time scale of a few femtoseconds (see [21 [3] and references 
therein). The projectile potential energy is given by the sum over the sequen- 
tial ionization energies up to the charge state q and can easily reach 5-50 keV 
for projectile charge states q ~ 20 - 50. Recently, the focus of research has 
shifted to applications in material science, where highly charged ions promise 
the possibility for well-controlled nano-sized surface modifications. Different 
topographies of surface modifications have been observed, depending on the 
projectile potential and kinetic energy and the target material [U |5] (see (6] 
for a recent review). 

Even for seemingly similar wide-bandgap insulator targets, namely ionic 
crystals of alkaline earth halides and alkali halides, experiments have found 
different types of surface modifications. On KBr, pits of one atomic layer 
depth have been observed [7] while on CaF 2 nano-meter sized hillocks pro- 
truding from the surface have been found [3]. In a recent contribution [lj, pits 
have also been found on CaF 2 surfaces if the sample was etched after irradi- 
ation ("etch pits"). The appearance of hillocks and etch pits as a function of 
both potential and kinetic energies of the impinging projectiles suggests the 
existence of three regions in parameter space: for low kinetic and potential 
energies, no surface modifications are observed while for a combination of 
larger kinetic and potential energies, only etch pits but no hillocks are found. 
Finally, for sufficiently large projectile potential energy both nano-hillocks 
and etch pits are observed only weakly dependent on the concomitant ki- 
netic energy. The appearance of these two different surface modifications for 
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the same target material was associated with the existence of two different 
thresholds, one for etchable damage and one for hillock formation. 

The involved complex many-body dynamics spanning many orders of 
magnitude in time pose a considerable challenge to theory. In this con- 
tribution, we present a first attempt of a multi-scale modelling making use 
of the different times scales rendering a qualitative, yet microscopic, model 
of surface modifications. 

The plan for this paper is as follows. We briefly review theoretical models 
for the transfer of projectile potential and kinetic energy to the electronic and 
subsequently to the lattice degree of freedom. Then, we present a proof-of- 
principle molecular dynamics simulation with the heated nanoscale surface 
region as initial condition that qualitatively reproduces the existence of two 
thresholds associated with the creation of etchable defects and surface re- 
structuring in the form of nano-hillocks in agreement to experiment. Atomic 
units are used throughout unless stated otherwise. 

2. Scenario for energy deposition 

The multi-scale model for the energy deposition by slow highly charged 
ions is based on the scenario in which the initial deposition into the the elec- 
tronic degree of freedom and subsequently into the lattice degree of freedom 
proceeds on shorter time scales than the subsequent atomic rearrangements 
and structural changes. This allows to treat such processes separately. 

As the slow highly charged ion approaches the surface, the primary energy 
conversion from potential energy of the projectile to electronic energy of the 
target proceeds by a complex de-excitation cascade, both above and below 
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the surface. As the ion approaches the surface, electron capture to highly 
excited projectile states sets in [8HTU]. followed by a subsequent de-excitation 
by Auger decay processes and radiative decay. Close to the surface, the self- 
image interaction leads to an upwards shift in energy of the highly excited 
states, leading to re-ionization of the projectile. As a consequence, a large 
fraction of projectiles are still highly charged upon entering the first few 
target layers. De-excitation below the surface continues with states below 
or in resonance with the target valence states, leading to moderately excited 
projectile states. They decay by Auger sequences emitting electrons with low 
to intermediate energies up to a few hundred eV. In the case of inner-shell 
holes (e.g., Ar 17+ ), a small number of keV electrons are released. The kinetic 
energy of the projectile determines the depth at which the neutralization 
sequence is complete, with penetration depths of a few nm for kinetic energies 
of < 10 5 eV jl]. The process taking place within the first ~ 20 femtoseconds 
of the interaction can be simulated using the classical-over-the-barrier model 
[8] and its extension to insulator surfaces [TT]. A quantitative analysis of 
slow highly charged ion impact on CaF 2 is given in [12], relating the charge 
state of the projectile to the electron emission spectrum around the impact 
area. 

These liberated electrons lead to the subsequent electron cascade, trans- 
ferring part of their kinetic energy to the target lattice by electron-phonon 
coupling. The relevant time scale for this process is a few hundred femtosec- 
onds. This electron cascade can be analyzed by a classical trajectory Monte 
Carlo simulation taking into account elastic and inelastic scattering events, 
leading to the creation of secondary electrons and to excitations of phonons 
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in the interaction with crystal atoms [HEI2]- One key result, somewhat coun- 
terintuitive at first sight, is that the low-energy part of the electron spectrum 
(E < 200 eV) is far more effective in the local heating of the target lattice 
than the energetic part (~ keV) of the electron spectrum. Even though their 
available energy is lower, their contribution to the heating of the lattice near 
the impact area is larger because they deposit their energy in close proximity 
to the impact site leading to a large energy density whereas high-energy 
electrons transport and spread their energy over a much larger volume due 
to their larger mean free path. One consequence of this rapid diffusion is 
that only a part of the projectile energy is effectively deposited as heat near 
the impact site. Typically, for an initial potential energy of the order of ~ 10 
keV (corresponding to Xe 24+ ) only 20 % are released as kinetic energy of 
low energy electrons, part of which is deposited as heat near the impact site. 
In total, an energy density near the melting point is reached in a region of 
diameter ~ 3 nm (fig. [T|. The influence of the kinetic energy on this process 
is primarily to control the shape and size of the heated "core" volume of the 
impact site. The shape varies from hemispherical for slow ions with smaller 
penetration depth to a flame-like shape for faster ions with larger penetra- 
tion depth. In addition and of less importance, the energy loss of the ions 
controlled by its nuclear stopping power contributes to the heating as well. 
At projectile energies of up to a few hundred keV (right panel in fig. [T|, the 
nuclear stopping power in the bulk is near its maximum of 2500 eV/nm. The 
energy is delivered via a collision cascade which is strongly directed towards 
the bulk. SRIM calculations [13] show that the bulk stopping power is only 
reached after several tens of nm inside the material where the deposited en- 
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ergy is spread over an area of several 100 nm 2 perpendicular to the projectile 
path leading to a low energy density. In the first few nm inside the solid, the 
local stopping power typically reaches only 20 % of the bulk stopping power. 
While the resulting energy density is too low to induce surface restructur- 
ing on its own it synergistically works together with the electron cascade to 
thermally activate the lattice. 

3. Simulation of target restructuring 

The physical scenario laid out in the preceding section results in strong 
local heating in a region near the impact point. We parametrize the outcome 
of the electron cascade by the linear size of the impact region, its radius R 
(assuming a hemisphere), and its effective temperature, T. These parameters 
serve as input for our proof-of-principle molecular dynamics (MD) simulation 
by which we explore the mechanisms leading to melting, etchability, and 
eventually surface restructuring. 

The target atoms are described as point particles interacting via static 
ionic pair potentials with short range admixtures of the Born-Mayer-Huggins 
form, 

Vij(r) = Z^/r + A^l™ - C tJ /r 6 (1) 

where we assume the bare ionic charges Zc a = +2 and Zp = — 1 using 
parameters from [13] which are tuned to give the correct lattice constant and 
F vacancy formation energy (Frenkel defect energy) at low temperatures. We 
neglect the preceding de-excitation and charge transfer processes and assume 
that all atoms are in their electronic ground state, giving and upper bound 
for crystal stability. 
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The size of the simulation box was up to 13.9x11.4x4.6 nm 3 (~ 55000 
atoms). Periodic boundary conditions are assumed perpendicular to the sur- 
face to avoid edge effects. The evaluation of Coulomb forces is performed in 
a shifted-force scheme ("Wolf method" [ToTflT] ) with a cut-off of 0.9 nm, the 
accuracy of which has been demonstrated for a variety of systems including 
ionic crystals [TTl [T8] . We use a Berendsen thermostat [19] to thermalize 
an ensemble of atoms. The outermost layers of the simulation box can be 
coupled to a heat bath at room temperature to mimic the influence of the re- 
maining crystal and act as a heat sink. The simulation reproduces a range of 
ground state and equilibrium properties of CaF 2 such as superionicity [20J, 
sound velocity, melting point, and thermal diffusivity [2TJ, which is a key 
quantity for the transport of heat away from the impact region. One ma- 
jor limitation of the present simulation is that the atomic potentials (eq. [T]) 
do not allow for charge transfer and transitions between potential hypersur- 
faces giving rise to chemical reactions. Therefore, the simulation results are 
intended to give primarily qualitative insights. 

The virtual crystal is first equilibrated at a temperature of 300 K for 3 
ps. To simulate the effect of a highly charged ion impact and the resulting 
heating of the target near the impact point, a central hemisphere (radius R) 
is heated to a temperature T within 100 fs, the time scale of the electron 
cascade (fig. [2]). A simulation snapshot after the heating is displayed in 
fig. [2] (b). We observe a large number of thermally activated F~. Within our 
simulation, most of the activated F ions find their way back to F vacancy sites 
as the system cools down (simulation time 15 ps, simulation snapshot in fig. [2] 
(c)), which is to be expected for a non-amorphizable material such as CaF 2 
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|22j . This reconstitution is, however, in part an artifact due to the absence of 
reactive scattering processes which are difficult to incorporate into a classical 
molecular dynamics simulation in a consistent way [23]. For example, the 
thermally activated F~ ions can form fluorine gas by the charge transfer 
reaction Ca 2+ + 2 — > Ca + F2. Similarly, neutral F° "holes" created in 
the F sub-lattice can covalently bond 2 F° — > F2, even without need for a 
nearby Ca. The newly formed F 2 molecules could aggregate as interstitials 
and evaporate, leaving behind Ca aggregates near the surface constituting 
the hillock, similar to bulk restructuring observed in CaF 2 after swift heavy 
ion impact [2H [25] . 

Depending on R and T, the formation of hillocks can be observed (fig. [2] 
(c)) as a patchy partial triple layer on the surface of the crystal. To quantify 
etchable defects we determine the number of thermally activated atoms of the 
target lattice by the number of F ions that end up further than half a nearest- 
neighbor distance of the F-sublattice (~ 0.14 nm) away from their original 
position since those atoms would have a significant chance to participate in 
chemical reactions. A large number of "activated" atoms weakens the target 
structure and leads to increased sensitivity to etching. 

We consider now the dependence of formation of hillocks and of etchable 
defects as a function of R and T describing the width and temperature of the 
heated zone. Fig. [3] shows contour maps of the simulated height of hillocks 
(panel (a), linear color scale) and the number of thermally activated atoms 
in the simulation box (panel (b), logarithmic color scale). For both observ- 
ables thresholds are identified (marked by black lines). Since both hillock 
height and the number of activated atoms are statistical quantities the re- 
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suits present an average over 10 simulation runs each on a grid of 5x5 points 
and interpolated in between. The two thresholds for hillock formation on one 
hand and etchable defects on the other hand are compared in fig. [3] (c). In 
the experiment, the hillock formation is always accompanied by etchability 
[1] , consistent with the present simulation where a large number of activated 
atoms is visible beneath a hillock (compare fig. [3^i, |3}d and right-hand panel 
offig.§. 

Typical snapshots after 15 ps for three different T at fixed R =3 nm 
(marked by points 1,2, and 3 in fig. [3]) show that no damage is visible for 
low temperatures (fig. [4] top panel), while a number of activated atoms, con- 
sisting mainly of F atoms, concentrated around the impact point appears for 
moderate temperatures (fig. Reenter panel). For high temperatures (fig. [4] 
bottom panel) the number of activated atoms is even higher now also includ- 
ing Ca atoms near the impact point. A hillock protruding from the surface 
is observed with a typical height of one triple layer ( ~0.3 nm). 

The two thresholds observed in the simulation are consistent with the ones 
observed in the experiment pj. For a more detailed quantitative connection 
the relationship between the experimentally controlled projectile parameters 
(potential and kinetic energies) with the MD simulation parameters R and T 
would be required. While the dependence of R and T on the potential energy 
can be approximately determined within our simulation for the initial heating 
by hot electrons, the dependence on the kinetic energy is more involved 
as increased energy dispersion (see fig. [T]) is partially offset by increased 
nuclear stopping. The combined effect of both the de-excitation cascade and 
the collision cascade is currently not included in our heat deposition model. 
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This exact relationship is, however, not of importance for the qualitative 
observation of two thresholds, which are crossed by any path leading from 
low values of (R,T) to large (R,T) fig. [3] (c). 

For future nano-structuring applications it is of interest to explore regions 
of larger T. The latter could, possibly, be reached by more highly charged 
ions (e.g. U 92+ ) and/or faster projectiles. With a five-fold increase in temper- 
ature (deposited thermal energy) compared with the one from nano-hillock 
formation discussed above such that the sublimation threshold is reached we 
observe a pronounced crater, i.e. a combined hillock and pit (fig. [5]). The 
depth of the crater is ~ 1.5 nm corresponding to a sputter yield of ~ 440 
atoms (volume of the crater). In addition, a swollen region forming the 
elevated crater rim (height ~ 1 nm) is visible. The total thermal energy 
deposited in the simulation is ~ 6 keV, which might be reached with the 
highest charge states or, more likely, with swift heavy ions. 

4. Conclusions and outlook 

In this contribution, we have theoretically investigated surface restruc- 
turing by highly charged ions. We have reviewed the mechanisms underlying 
the surface restructuring and the theoretical tools to model the deposition of 
the potential energy of the highly charged ion to the target lattice. Within 
a proof-of-principle molecular dynamics simulation to estimate the effects of 
the HCI impact on the target lattice we could identify the appearance of 
two thresholds associated with etchability and with nano-hillock formation 
in qualitative agreement with experiment. 

To improve the current theoretical model, the inclusion of charge transfer 
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reaction processes and chemical changes during the de-excitation sequence 
is of paramount importance. While the formation of surface modification 
happens on a time scale of a few picoseconds, the chemical structure of the 
surface and bulk continues to change through defect agglomeration and defect 
mediated desorption on a longer time scale. 

Along these lines, measurements of angle- and species-resolved sputter 
yields could provide important clues as to the microscopic dynamics of the 
nano-melt during hillock formation. An in-situ measurement of the chemical 
composition of hillocks would give information about the final state of the 
system as well as a more complete characterization of hillock formation for 
potential applications in nano-science and could provide important bench- 
marks for MD simulations that include reactive processes. 
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Figure 1: Simulated thermal energy deposition due to impact of Xe 33+ for low (left) and 
larger (right) projectile kinetic energy. For details see text and [I]. 
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Figure 2: Simulation of etchable defects and hillock formation, (a) Shape of the hemi- 
spherical region (radius R) heated to a temperature T in a CaF 2 crystal (grey: Ca, yellow: 
F). (b) Snapshot after the heating (duration 100 fs). (c) The cooled crystal after 15 ps 
simulation time shows few dislocations in the bulk and a patchy triple layer on the surface 
is visible as a hillock under the AFM. 
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Figure 3: Simulation results as a function of MD parameters R and T respectively describ- 
ing the size and temperature of the thermally activated region created by the preceding 
HCI impact, (a) Height of simulated hillocks, (b) Number of thermally activated atoms 
surmised to lead to etchable defects (details see text), (c) Combined graph showing two 
different thresholds, for example along the R = 3 nm contour marked by points 1, 2 and 
3 corresponding to simulation results shown in fig. [4j 
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Figure 4: Typical simulation results showing the transition from no restructuring (top) to 
a large number of activated atoms near the impact region (center) to surface modification 
(bottom). Light atoms are F while darker atoms are Ca. The distance of each atom to 
its original position at the beginning of the simulation is overlayed as color code. The 
simulations correspond to R = 3 nm and T = 1000 K, 2000 K and 3000 K respectively, 
marked by point labels 1, 2 and 3 in fig. |3jc). Figure created by AtomEye [26] . 



18 



Figure 5: Snapshot of MD simulation displaying crater formation. The maximum tem- 
perature of the central hemisphere (radius R = 2 nm) was taken as T = 18000 K. Left: 
Top-down view, right: cut through center, with color according to the distance from the 
surface layer. System size is 14.7x13.4x4.5 nm 3 (68400 atoms). Figure created with 
AtomEye [26]. 
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